log using "$Rep_smokelabor/2_analysis/output_logs/appendix_figure3.log", replace
	
	use "$Rep_smokelabor/1_build/regdata/county_day.dta", clear
	
	tsset countyfip rfrnc_dt
	
	foreach poll in pm25 o3 so2 co no2 pm10 {
		reghdfe std_`poll' F(20(-1)1).hms_deep hms_deep L(1(1)20).hms_deep, a(county_doy st_yr) cluster(COUNTY10 rfrnc_dt) 
	}
	
	** plot
	import excel using "$Rep_smokelabor/2_analysis/output_figures/appendix_figure3.xlsx", clear sheet("std") first
	
	tw 	connected pm25 eday, col(blue) msize(0.6)  ///
		xtitle("Days since smoke exposure") ytitle("Change in concentration (S.D.)") ///
		xline(0, lcol(black) lp(dot) lw(0.6)) ///
		yline(0, lcol(black) lp(dot) lw(0.6)) ///
		ylab(-0.1(0.1)0.4, nogrid) ///
		text(0.4 20 "PM2.5", size(vlarge) place(sw) box just(left) fcolor(white) lw(0.4) margin(medsmall)) ///
		legend(off) ///
		graphregion(col(white))
	gr export "$Rep_smokelabor/2_analysis/output_figures/appendix_figure3_pm25.pdf", replace
	
	tw 	connected pm10 eday, col(blue) msize(0.6) ///
		xtitle("Days since smoke exposure") ytitle("Change in concentration (S.D.)") ///
		xline(0, lcol(black) lp(dot) lw(0.6)) ///
		yline(0, lcol(black) lp(dot) lw(0.6)) ///
		ylab(-0.1(0.1)0.4, nogrid) ///
		text(0.4 20 "PM10", size(vlarge) place(sw) box just(left) fcolor(white) lw(0.4) margin(medsmall)) ///
		legend(off) ///
		graphregion(col(white))
	gr export "$Rep_smokelabor/2_analysis/output_figures/appendix_figure3_pm10.pdf", replace
	
	tw 	connected o3 eday, col(blue) msize(0.6)  ///
		xtitle("Days since smoke exposure") ytitle("Change in concentration (S.D.)") ///
		xline(0, lcol(black) lp(dot) lw(0.6)) ///
		yline(0, lcol(black) lp(dot) lw(0.6)) ///
		ylab(-0.1(0.1)0.4, nogrid) ///
		text(0.4 20 "O3", size(vlarge) place(sw) box just(left) fcolor(white) lw(0.4) margin(medsmall)) ///
		legend(off) ///
		graphregion(col(white))
	gr export "$Rep_smokelabor/2_analysis/output_figures/appendix_figure3_o3.pdf", replace

	tw 	connected co eday, col(blue) msize(0.6) ///
		xtitle("Days since smoke exposure") ytitle("Change in concentration (S.D.)") ///
		xline(0, lcol(black) lp(dot) lw(0.6)) ///
		yline(0, lcol(black) lp(dot) lw(0.6)) ///
		ylab(-0.1(0.1)0.4, nogrid) ///
		text(0.4 20 "CO", size(vlarge) place(sw) box just(left) fcolor(white) lw(0.4) margin(medsmall)) ///
		legend(off) ///
		graphregion(col(white))
	gr export "$Rep_smokelabor/2_analysis/output_figures/appendix_figure3_co.pdf", replace
	
	tw 	connected no2 eday, col(blue) msize(0.6)  ///
		xtitle("Days since smoke exposure") ytitle("Change in concentration (S.D.)") ///
		xline(0, lcol(black) lp(dot) lw(0.6)) ///
		yline(0, lcol(black) lp(dot) lw(0.6)) ///
		ylab(-0.1(0.1)0.4, nogrid) ///
		text(0.4 20 "NO2", size(vlarge) place(sw) box just(left) fcolor(white) lw(0.4) margin(medsmall)) ///
		legend(off) ///
		graphregion(col(white))
	gr export "$Rep_smokelabor/2_analysis/output_figures/appendix_figure3_no2.pdf", replace
	
	tw 	connected so2 eday, col(blue) msize(0.6) ///
		xtitle("Days since smoke exposure") ytitle("Change in concentration (S.D.)") ///
		xline(0, lcol(black) lp(dot) lw(0.6)) ///
		yline(0, lcol(black) lp(dot) lw(0.6)) ///
		ylab(-0.1(0.1)0.4, nogrid) ///
		text(0.4 20 "SO2", size(vlarge) place(sw) box just(left) fcolor(white) lw(0.4) margin(medsmall)) ///
		legend(off) ///
		graphregion(col(white))
	gr export "$Rep_smokelabor/2_analysis/output_figures/appendix_figure3_so2.pdf", replace

log close